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Abstract 

We directly sample the exponential moments of the canonical phase for various 
quantum states from the homodyne output. The method enables us to study the 
phase properties experimentally, without making the detour via reconstructing 
the density matrix or the Wigner function and calculating the phase statistics 
from them. In particular, combing the measurement with a measurement of the 
photon-number variance, we verify fundamental number-phase uncertainty. 
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Optical homodyne tomography has been a powerful method for quantum- 
state measurement, because the measured quadrature-component distributions 
p(x, $) = (x, fl\g\x, $) contain all the information about the quantum state [\x, i?) 
being the eigenstates of = (e _M? d + e^a))]. In this way, the Wigner 
function Q and the density matrix H |], [| can be inferred from experimental 
data. When the state (in a chosen representation) is known, then the mean values 
of arbitrary quantities can be calculated. 

Each measurement is associated with errors (at least, the statistical error) 
which propagate in the calculation procedure. Hence, the error of a calculated 
quantity can be too large to be acceptable. For example, if we are interested in 
photon-number moments (n k ), we have to multiply the diagonal density- matrix 
elements by n h and perform the sum. Since the statistical error of g n>n does 
not vanish with increasing n, the total error is infinite after summing all the 
terms. Truncation of the sum can avoid this trouble for sufficiently low-order 
moments, whereas for high-order moments also small values of Q n ^ n with large n 
may be essential. To overcome this problem, it has been suggested to sample the 
desired quantities directly from the measured data, without the detour via the 
density matrix (or other state representations). Formulas have been derived that 
are suited for direct sampling of normally ordered moments of photon creation 
and destruction operators ||, and an extension to arbitrary quantities that admit 
normal-order expansion has been given ||. Quite recently sampling functions for 
measuring the exponential moments of the canonical phase have been derived 0. 

In this paper we apply the sampling method to an experimental determina- 
tion of the exponential moments of the canonical phase for various single-mode 
quantum states. The results are then used for determining the phase distribution 
as Fourier transform of the exponential phase moments. In particular, the sam- 
pled first exponential phase moment already defines a phase uncertainty || ||. 
Using the (simultaneously) sampled photon-number variance a verification of the 
number-phase uncertainty predicted theoretically in |9j is given. 

The exponential phase moments are defined by the Fourier components of 
the phase distribution P(ip), i.e., ^k = 1^ dtp e lk(p P{ip>). For the canonical phase 
they are given by fT(| ^>k = (E k ) if k > 0, and = *- if A; < 0, where E = (n 



+ 1) 1 ' 2 a, h = a) a being the photon-number operator. It can be shown that 
can be obtained from the quadrature-component distributions p(x, #) as |7] 

/ d# dxK k (x,$)p(x^), (1) 

where Kk(x,d) is a well-behaved integral kernel suited for direct sampling of ^ 
from the homodyne output for any normalizable quantum state. 

Knowing ty^, the phase distribution P(<p) can be obtained according to P(<p) 
= (27r) _1 Y^"kL-oo €~ lkip ^k- However, the first moment already contains essential 
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information about the phase properties. It can be used to introduce a mean 
phase ip = arg^x and a phase uncertainty Aip = arccos|\I/i| , which implies a 
number-phase uncertainty relation || 

An tan Ay? > |. (2) 

Note that for the number- uncertainty An=((n 2 ) — (n) 2 ) 1 / 2 the quantities (h) and 
(n 2 ) can also be obtained by direct sampling according to a relation of the form 
(HD, with the integral kernel being given in ||. Hence, homodyne detection can 
be regarded as the most direct way that has been known so far for experimental 
verification of the uncertainty relation fl2|). 

The experimental setup is the same as in |4j]. Its central unit is a monolithic 
standing-wave lithium-niobate optical parametric amplifier (OPA) pumped by a 
frequency-doubled continuous- wave Nd:YAG laser (1064 nm). Operated below 
threshold, the OPA is a source of squeezed vacuum. We study the spectral 
components of the field around a frequency offset by Q/2tt= 1.5 or 2.5 MHz from 
the optical frequency u, to avoid low-frequency laser excess noise. To generate 
bright light, a very weak wave split off the main laser beam is phase-modulated 
by an electro-optic modulator (EOM) at the frequency fl and injected into the 
OPA through its high reflector port. The carrier frequency us is kept on- resonance 
with the cavity and the two "bright" sidebands uj ± Q are well within the cavity 
bandwidth. By turning the modulation off, we obtain squeezed vacuum, by 
blocking the OPA pump wave, we are left with coherent excitations. 

The signal is analyzed at a homodyne detector, whose output current is mixed 
with an electrical local oscillator phase-locked to the modulation frequency, low- 
pass filtered and recorded with a high speed A/D converter. Since the squeezed 
states are essentially two-mode states, a two-mode detection is crucial for ob- 
taining the correct statistics of the light field. We remark however, that this 
type of measurement may need modifications for general states of the light field. 
The quadrature-component distributions p(x, i?) are obtained by subdividing the 
recorded noise traces into 128 equal length intervals and subsequently forming 
histograms of 256 amplitude bins, normalizing the absolute bin width using as 
reference the distribution of a vacuum state. 

In Figs. [T]-^, the sampled exponential phase moments k — 1, 2, . . . , 20, 
are shown for a phase-squeezed state (Fig. [I]), a state squeezed at a phase angle 
of 48° (the difference between the argument of the squeezing parameter and the 
argument of the displacement parameter) and a squeezed vacuum (Fig. |3|). The 
error bars indicate the statistical error. Since the main source of inaccuracy is the 
fluctuation of the local oscillator, the error is dominated by the systematic one. 
The canonical phase distributions obtained from the sampled moments (Fourier 
components) are plotted in Fig. |j. Since in Fig. [I] |^| decreases with increasing 
|fc| slower than in Fig. [l|, the phase distribution of the phase-squeezed state in 
Fig. |](a) is more sharply peaked than that in Fig. f|(b). Figure. |](c) clearly 
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reveals the double-peak structure of the phase distribution of a squeezed vacuum. 
Note that the small oscillations in the figures (which also include negative values) 
mainly result from systematic errors. 

Examples of (p and A<p together with the measured mean photon number 
n = (h) and photon-number uncertainty An are given in Tab. |1] for various 
states prepared in the experiment. The last row shows the resulting values of 
the number-phase uncertainty product AntanA<£>, which are in full agreement 
with the predicted inequality The (near-)coherent states (A,B) and phase- 
squeezed states (E,F) are seen to exhibit relatively small phase uncertainties. 
Note that the coherent state (B) has the smallest phase uncertainty. Relatively 
large phase uncertainties are observed for the amplitude-squeezed states (C,D) 
and the state (H) squeezed at a phase angle of 48°. As expected, the near- 
maximum phase uncertainty A(pmir/2 corresponds to the squeezed vacuum (G) 
(the "ellipse" in the phase space is centred at the origin). Therefore, for this 
state the uncertainty product An tanAyj achieves a very large value. The small- 
est value of the uncertainty product is observed for the coherent state (A). It is 
close to its limit 1/2. With respect to the photon number, the amplitude-squeezed 
state (C) is seen to be sub-Poissonian. 

In summary, we have sampled the exponential moments of the canonical phase 
directly from the homodyne output for various coherent and squeezed states 
produced in a continuous optical field by means of parametric amplification. This 
has enabled us to study the canonical phase statistics experimentally, without the 
necessity of state reconstruction, which saves calculation effort and reduces the 
statistical error. In particular, from the sampled first-order exponential phase 
moment and the simultaneously sampled first- and second order photon-number 
moments we have determined phase and number uncertainties and shown that 
the uncertainty products are in agreement with the theoretical prediction. 

This work was supported by the Deutsche Forschungsgemeinschaft. 
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state 


A 


B 


C 


D 


E 


F 


G 


H 


<P 


0.02 


1.59 


3.13 


3.13 


-3.13 


1.56 


2.62 


0.78 


Aip 


0.17 


0.14 


0.30 


0.54 


0.17 


0.14 


1.56 


0.31 


n 


8.62 


25.94 


9.19 


14.72 


8.95 


38.45 


6.92 


26.46 


An 


3.03 


6.20 


2.30 


7.21 


5.38 


25.05 


10.02 


14.85 


An tanAy? 


0.52 


0.87 


0.71 


4.32 


0.92 


3.53 




4.75 



Table 1: Measured values of <p, Aip, n= (n), and An, and the resulting values 
of the number-phase uncertainty product An tan Aip for various quantum states 
[(A,B) coherent states; (C,D) amplitude-squeezed states; (E,F), phase-squeezed 
states; (G) squeezed vacuum; (H) state squeezed at a phase angle of 48°]. 
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Figure 1: Measured exponential phase moments for a phase squeezed state 
[state (F) in Tab. |TJ [(a) real part; (b) imaginary part]. 
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Figure 2: Measured exponential phase moments for a state squeezed at a 
phase angle of 48° [state (H) in Tab. |I|] [(a) real part; (b) imaginary part]. 
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Figure 3: Measured exponential phase moments for a squeezed vacuum [state 
(G) in Tab. [l| [(a) real part; (b) imaginary part]. 
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Figure 4: The canonical phase distribution P(<p) reconstructed from 20 measured 
exponential phase moments given in Figs. [T| - |3] is shown for (a) a phase- 
squeezed state [state (F) in Tab. [L|, (b) for a state squeezed at a phase angle of 
48° [state (H) in Tab. U and (c) for a squeezed vacuum [state (G) in Tab. |]. 



